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O ! Abstract 



The structure of a relativistically hot, strongly magnetized jet is investigated at large distances 
from the source. Asymptotic equations are derived describing coUimation and acceleration of the 
externally confined jet. Conditions are found for the transformation of the thermal energy into the 
fluid kinetic energy or into the Poynting flux. Simple scalings are presented for the jet collimation 
angle and Lorentz factors. 
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^ ; 1 Introduction 

Highly collimated, relativistic outflows are commonly observed in compact astronomical objects 
> ■ 

I (gamma-ray bursts, active galactic nuclei, galactic black- hole and neutron-star binaries). It is 
l/^ . widely believed that in all these cases, the jet is driven by rotating, twisted magnetic fields (see, 
^ ' e.g., recent reviews in ref.pj). The rapidly spinning central body (neutron star, accretion disk, 
black hole) twists up the magnetic field into a toroidal component and the plasma is ejected by 
the magnetic tension. In relativistic flows, the energy per particle significantly exceeds the rest 
mass energy therefore in order to create a relativistic jet, one has to convey a significant energy 
^ ' to a small amount of the matter. The main advantage of the magnetic launch mechanism is that 
- - . the magnetic field lines, like driving belts, could in principle transfer the rotational energy to a 
low density periphery of the central engine thus forming a baryon pure but energetic outflow. The 
relativistic velocities could potentially be achieved if the magnetic energy density in the plasma 
frame exceeds the plasma energy density. In such outflows, the energy is transported, at least 
initially, in the form of the Poynting flux. The question is how and where the electromagnetic 
energy is eventually transformed to the plasma energy. 

In the scope of ideal MHD, the energy could be transferred to the plasma only via gradual 
acceleration by electromagnetic stresses. In non-relativistic flows, the plasma is accelerated cen- 
trifugally when sliding along the rotating poloidal field lines. The azimuthal field is generated 
only when the plasma inertia becomes comparable with the magnetic stresses so that the field 
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lines have to bend backwards. This imphes that a good fraction of the Poynting flux is converted 
into the kinetic energy of the flow already when the azimuthal field becomes comparable with the 
poloidal one. In relativistic, highly magnetized flows, the magnetic force is generally balanced not 
by inertia but by the electric force. The azimuthal field becomes comparable with the poloidal one, 
both of them being comparable with the electric field, at the light cylinder defined as a surface 
on which the corotational velocity is equal to the speed of light (in differentially rotating magne- 
tospheres, this surface is not a cylinder but we retain the standard term, which has come from 
the pulsar theory). The fluid kinetic energy remains small at the light cylinder. Beyond the light 
cylinder, the conservation of the magnetic flux implies that the poloidal magnetic field decreases 
as where r is the cylindrical radius of the jet. The azimuthal field and the electric field 

decrease only as 1/r being close to each other so that the fluid is only slowly accelerated by a small 
residual force. By this reason, the acceleration zone is extended well beyond the light cylinder 
so that formation of relativistic jets spans a very large range of scales. In recent studies, both 
numerical and analytical [21 [Sj HI [5l |6l [TJ [8] , the general conditions were formulated for coUimation 
and acceleration of relativistic MHD jets and the efficiency of the Poynting flux into the kinetic 
energy conversion was thoroughly examined. Non-steady Poynting dominated outflows have also 
being studied [21 HDl IHl 112] • In these works, only cold flows have been addressed. Here we relax 
this assumption and study relativistically hot, highly magnetized jets. By relativistically hot we 
mean the fluid with the pressure exceeding the rest mass energy density, p > pc^. Such outflows 
are believed to be formed in gamma-ray bursts (see, e.g., reviews [T3| [H]). Here we study such 
outflows in the far zone where most of acceleration occurs. 

When the flow expands, the fluid cools down. In the pure hydrodynamical case, the thermal 
energy is converted into the kinetic energy of the fluid. It follows immediately from the Bernoulli 
equation that the relativistically hot fluid with the adiabatic index F = 4/3 is accelerated such that 
the flow Lorentz factor grows proportionally to the jet radius. In a magnetized flow, the thermal 
energy could be converted not only into the kinetic energy but also into the Poynting flux. The aim 
of this study is to explore the fate of the thermal energy in expanding, Poynting dominated jets. 
We consider outflows confined by the pressure of the external medium because only in this case 
the jet could be collimated. The paper is organized as follows. In section 2, the general equations 
governing ideal relativistic axisymmetric flows are presented. In section 3, asymptotic theory is 
developed describing such outflows in the far zone. In sect. 4, the theory is applied to narrow jets 
from rigidly rotating sources. Conclusions are presented in sect. 5. 
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2 Relativistic MHD equations for a steady axisymmetric 
flow 



For the sake of consistency and in order to introduce notations, let us shortly review the basic 
theory of relativistic, magnetized outflows p3| [T6| [T7] . In order to describe ideal MHD flows one 
has to employ the Maxwell equations, the flux-freezing condition and the continuity, entropy and 
momentum equations. For relativistic, steady state flows these equations are written as 

V-B = 0; VxE = 0; (1) 

V X B = 47rj; V ■ E = A-npe] (2) 
E + V X B = 0. (3) 

V ■ (p7v) =0; V ■ V (^-^^ = 0; (4) 
P7(v ■ V)(/i7v) = -Vp + PeE + j X B. (5) 

Here v is the velocity of the flow and 7 = (1 — f ^)~^/^ the associated Lorentz factor, E and B the 
electric and magnetic fields as measured in the lab frame, pe and j the charge and current densities 
also in the lab frame, and p, h and p the plasma proper mass density, specific enthalpy and the 
pressure. The speed of light is taken to be unity throughout the paper. The above equations 
should be complemented by the equation of state, which we assume to be polytropic with the 
index F. Then the specific enthalpy may be presented as 

.^1 + ^^. (6) 

Here the first term represents the rest mass energy. 

In axisymmetrical flows, the magnetic field is conveniently decomposed into the poloidal and 
azimuthal components. Bp and -B<^. The poloidal field is expressed via the flux function as 

Bp = ^V^ X 0. (7) 

We use cylindrical (r, 0, z) coordinates; the hat denotes unit vectors. The flux freezing condition 
^ implies that the magnetic surfaces are equipotentials so that one can write 

E = rn{^)Bp; (8) 

where fl is the angular velocity of the magnetic field line. The position of the light cylinder is 
found from the condition Qr = 1. The plasma streams along the magnetic surfaces so that one 
can decompose the fluid velocity into the poloidal (along Bp) and azimuthal components. Then 
Eqs. ([3D and ([8]) yield 

BpV^ - B^Vp = rfi(^')5p; (9) 
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which imphes that the fluid shdes along the magnetic field lines rigidly rotating with the angular 
velocity Q. It follows immediately from Eqs. ([8]) and (Q that well beyond the light cylinder, 
Qr ^ 1, the electric and azimuthal magnetic field are close to each other and both are much larger 
than the poloidal field, E QrBp. 

If the shape of the flux surfaces is given, the flow is completely described by five integrals of 
motion. The first one is the angular velocity of the field line, fi(\E')), defined by Eq. Three 
others are the mass, energy and angular momentum-to-magnetic flux ratios, 

vm = (10) 

-Dp 

m = /.7 - ^ (11) 

V 

im = h^rvs-'^. (12) 
V 

The fifth integral is the specific entropy 

m = J- (13) 

These equations are obtained by projecting Eqs. (jlj) and (ISl) onto the direction of the flow, 
1 = n X 0, where n = V\E'/|V\E'| is the normal to the flux surfaces. The integrals of motion are 
generally found by prescribing the conditions at the inlet of the flow and from the condition of 
smooth passage of the solution through the critical points of the equations [171 [HI [HI [20] . Here 
we consider the flow in the far zone, outside of all critical points, therefore all integrals of motion 
are considered in this paper as given functions of \E'. 

Making use of the above conservation laws, one can express any physical quantity via, e.g., 7, 
h and ^ . For example, 

^0 = 7^ 1 - ■ (14) 



If the rotation velocity at the origin of the outflow is well below the speed of light, firin, f(^,in ^ 1, 
one can write Eq. iHM in the form 

-^ = 7i:fi-%l; (15) 



f2r \ h'j 

where the index "in" is referred to the parameters of the injected plasma. Below we will use this 
formula because it results in a bit simpler expressions than the general formula (IT^ . In the general 
case, we should just substitute /iin7in hj fi — Ql in all the following expressions. One can also found 
straightforwardly Vp as a function of 7, h and ^. Then the identity f ^ + f ^ + 7^^ = 1 may be 
written as the Bernoulli equation 



p 



2 



^ 1 /lin7i, 



1 2 



^2^2 ^ h'~) J 7' 



/ \ 2 

1 hmlni\ . 1 



1-^ +_ = 1; (16) 



which connects the Lorentz factor of the flow with the geometry of the flux tube defined by the 
flux function 

In order to get a complete set of equations, one has to make projection of the equation of 
motion onto the normal to the flux surface. This yields the transfield force-balance equation (the 
generalized Grad-Shafranov equation) 



1 

n 



n ■ r 



— n ■ Vp 



:n ■ V [r2(5' - E^' 



(17) 



Here 7^ is the local curvature radius of the poloidal field line (defined such that 7^ is positive when 
the flux surface is concave so that the coUimation angle decreases) 



_ ^ _n.(l ■ V)l = fi-[lx ((V X 1))] = -0-(V X 1). (18) 

The transfield equation ( !T7|1 is reduced to an equation for 7, h and \E' with the aid of Eqs. (|6]) and 
(ED - ([H. 



3 Relativist ic jet in the far zone 

In this section, we derive asymptotic equations describing the flow at large distances from the 
source. For cold relativistic flows the procedure was developed in ref. here we present the 
straightforward generalization to the hot case. We are going to obtain the Bernoulli and transfield 
equations in the limit fir ^ 1, 7 ^ 1. These two conditions are not independent for outflows that 
are initially Poynting dominated because such outflows are generally accelerated therefore such 
outflows are in any case highly relativistic in the far zone. 

The zeroth approximation in (f2r)~^ and 7"^ to the Bernoulli equation is obtained by taking 
the limit fir, 7 —t- 00 in Eq. f|T6|) . which yields 

ry(/i - h-i) = n^r^Bp. (19) 

This equation could be written, with account of Eqs. and f lTT]) . as 

B^ + E = 0. (20) 

In the transfield equation f lT7|) . the leading order in 1/r and I/7 terms are those in the right-hand 
side, n- V(i?0 — i?^), because the terms in the left-hand side are small either as Bp/B^ ~ {Qr)~^ or 
as r/TZ (the curvature of the flow lines in the far zone is small). This means that if one substitutes 
Eq. (120|) into the right-hand side of the transfield equation, one would kill the leading order terms. 
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The correct procedure [211 [7] is to expand the BernouUi equation (|T6l) to the first non-vanishing 
order in 1 /r and 1 /7 and only then to eliminate the leading order terms from Eq. ([1] 
Expanding Eq. f|T6l) yields 

„ fiV52 r /A. \2 



r 



1 



^inTii 



(21) 



Substituting this relation into the right-hand side of Eq. ( IT71) . one comes to an equation that does 
not contain terms which nearly cancel each other. Therefore one can safely simplify this equation 
further out. Namely, we can now neglect Bp as compared with substitute unity instead of Vp 
and also use the Bernoulli equation in the form f|T9|) in order to transform terms. The procedure is 
the same as that was described in sect. 4.1 of ref. |7]. The resulting asymptotic transfield equation 
is written as 



ri^iBp ^ rjBp 



n 



2r2 



1,2 2 

fi - 2/ii,7i„ + ^ 
7/1 



n ■ r 



n ■ Vp 



n ■ V 



1^272 



1 + 



^in7in 

hVtr 



(22) 



In the most interesting case of coUimated flows, 2; ^ r, one can take n - r = r and n - V = d/dr. 
When looking for the shape of the magnetic surfaces, one can conveniently use the unknown 
function r(\E',2) instead of \E'(r, 2;). Then, e.g., 

1(9^ 



Br, 



-IV^I 



r dr 



dr \ 



-1 



(23) 



In the same approximation, the curvature radius may be presented as (note that TZ is defined to 
be positive for concave surfaces) 

(24) 



1 

n 



Now the transfield equation for the coUimated flows in the far zone could be written as 




+ 



2/lin7 



h? '7-2 

in _|_ "■in /m 



dp 



hQr 



(25) 



This equation should be supplemented by the zeroth order Bernoulli equation (|T9l) . which is written 
in the new variables as 

(26) 



, ,dr ^]2 
(/i- /i7)— = —r- 



In order to get a closed set of equations, one has to express the pressure and the enthalpy via 
r and 7. In the far zone, Vp 1, <^ 1, the continuity equation (jl]), together with Eqs. ([9]) and 
(HH), gives 

rjBp _ ( r]_\'^ n-hj 



P 



47r7 



6 



47r7 



(27) 



Now one can write the entropy equation (1131) as 



jj, - hia7_ 

whereas the equation of state ([6]) is written in the form 

h-1 f r2^7in /i - /i7 



(29) 



/lin - 1 V ^^7 At - /lin7ii 

Equations (l25l) and f l26l) . supplemented by Eqs. ( l28l) and ( l29l) . form a complete set of equations 
describing relativistic, hot, magnetized, coUimated outflows in the far zone. We believe that these 
equations suit well to numerical solution because they do not contain terms that nearly cancel each 
other. However they also could be used in order to find simple analytical scalings in all limiting 
cases. 

Here we are interested in outflows subtending a finite magnetic flux \&o- If the flow is confined 
by the pressure of the external medium, pey±{,z)^ the pressure balance condition should be satisfied 
at the boundary: 

^{B^-E')+P = pU^). (30) 
ovr 

Making use of Eqs. f[T^ and fl^T]) . one writes the boundary condition in the form 

r/2(;i - h-i) 



1 



hVLr 



+ P) =Pey,t{z). (31) 



4 Relativist ically hot Poynting dominated jet 

In this section, we consider in detail a relativistically hot flow, 

/i > 1. (32) 

which was launched as Poynting dominated, 

hnlin < (33) 



When the flow expands, the plasma cools so that eventually the condition fl32|) is violated. The 
condition fl33|) ensures that the flow still remains Poynting dominated at this point therefore the 
condition 

h-i<^H (34) 
is fulfilled everywhere in the domain of interest. The enthalpy equation f l29p is now reduced to 



(35) 



From here on we will use the relativistic adiabatic index F = 4/3. The expanding accelerating 
fluid cools down so that that the condition (l32l) could be written as 

« /^fn^:^7m. (36) 

For simplicity, we consider the jet with a constant angular velocity, Q = const and homogeneous 
injection, rj = const; 7in = const. In this case, one can conveniently use the dimensionless variables 

X = fir; Z = Vtz. (37) 

We also assume that the energy integral, /i(^l/), is described by a linear function 

,(*)^/..,,„(l + |); (38) 

The last represents in fact the expansion of to the first order in \E' so that this expression is 
universally correct close to the axis of the fiow (see sect. 3 in ref. [7]). Numerical simulations show 
that in jets with a constant angular velocity, the linear function fl38p is a good approximation to 
/i all the way to the boundary of the fiow \E' = \E'o (see, e.g., fig. 4 in ref. |2j). In this expression, 
the first and the second terms describe the kinetic and the Poynting energy fiux, correspondingly. 
Note that the Poynting fiux goes to zero at the axis therefore the fiow is Poynting dominated only 
at ^ \E'; in this domain, the energy integral is reduced to 

. , 2fi2^ , , 

//(vj/) = . (39) 

4.1 The governing equation 

We are going to solve Eqs. ( |25|) and (!26|) under the conditions (132|) and The condition ( IM|1 
implies that one can neglect the second term in curly brackets in the left-hand side of the transfield 
equation (l25l) . The third term in the curly brackets is small at the inlet of the fiow however, the 
product h'-y could in principle grow with the distance. By virtue of Eq. (I35il this is possible only 
if 7 grows faster than r. We will see that this never happens therefore we neglect this term. 

Let us also demonstrate that the last term in the left-hand side of Eq. fl25|) is small as compared 
with the right-hand side of this equation. Making use of Eq. f[T^ one can estimate the ratio of 
these terms as P7^/ {VL^r'^B'^). The expression in the numerator is of the order of the plasma energy 
fiux whereas the expression in the denominator is, by virtue of Eqs. ([8]) and (!20|) . of the order of 
the Poynting fiux. Therefore the above ratio is small. 

Now the transfield equation is written, with the aid of Eq. (15^ . as 

„2 \ 2/3" 




1 + 



77i, 



(40) 



In the same approximation, the boundary condition ( pTj) is reduced to 



1 + 



77in 

XX2 



2/3' 



8vrpext(^)- 



(41) 



'I' = *0 



At the condition (IMI) and with account of Eq. (15^ one reduces the BernouUi equation 
an equation for X(^, Z): 

The general solution to this equation may be presented as 



to 



(42) 



X 



(43) 



where R{z) is an arbitrary function, which is in fact the radius of the jet. One sees that the 
structure of collimated, Poynting dominated jets is generally self-similar. Recall that this equation 
is valid only at \E' ^ \E'. Note that Eq. (l43l) implies that the poloidal magnetic field is homogeneous, 
oc X^; dBp/dX = 0; such a distribution is indeed seen in numerical simulations [1]. 
In order to find the function R{z), let us substitute Eq. (H3!) into the left-hand side of Eq. (HOjl 
and integrate the obtained equation between \1/ and \l/o: 




3 dZ^ 



1-] 

Ry 



(44) 



77in 



2/3' 



One can take the expression in the right hand side only at the upper limit because we will see 
that 7 grows not faster than X oc a/1' so that this expression grows with Making use of the 
boundary condition f HT]) one gets a closed equation for the jet shape: 



1 d'^R 



-3« 



dZ^ ^ i?2 



STipe^tR^ 



Presenting the external pressure as 



(45) 



(46) 



where po = Pext(l) is the external pressure extrapolated to the light cylinder one presents this 
equation in a dimensionless form 



We introduced the parameter 



fR 
dZ^ 



(3 



— +pViZ)R = 0. 



R2 ' 
^0 



(47) 



(48) 
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where Bq = Q^'^o is the characteristic magnetic field at the hght cyhnder. 

Having found the jet shape, one can easily find the full structure of the jet therefore Eq. fj47|) 
may be called the governing equation. If R{Z) is found from the governing equation, the shape 
of all the fiux surfaces is immediately obtained from Eq. fH3|) . In order to find the Lorentz factor 
of the fiow, one can substitute Eq. ( H3l) into Eq. (HOl) and then perform integration from to 
which yields an equation for 7 

It is interesting that the governing equation (H7|) coincides, to within change of variables, with 
that for cold fiows [3] 17]. Therefore the results for cold jets could be straightforwardly generalized 
to the hot case. If the external pressure is distributed according to a power law, 

^=^' (50) 

the general solution to the governing equation f H7|) could be found analytically (see ref. [7]). Below 
we describe these solutions. Since behavior of the solutions at large Z depends on the sign of k — 2, 
we consider different cases separately. 



2/3 



(49) 



4.2 The case k < 2. 

If K < 2, the fiow expands asjU El Ej 



3^«xi/4 



^=(^) ■ (51) 



This asymptotics could be found directly from Eq. ( 1471) by neglecting the first term. Substituting 
Eq. ( 15T|) into Eq. p7|) . one sees that at k < 2, the first term is really less than the second 
one. According to Eq. fHSj) . all fiux surfaces are similar to the surface of the jet boundary 
therefore in the transfield equation f l25l) . one can also neglect the term d'^r/dz^. Then the transfield 
equation becomes an ordinary differential equation describing cylindrical equilibria. This means 
that the structure of the jet at any z is the same as the structure of an appropriate cylindrical 
equilibrium configuration, parameters of the configuration being adiabatically changed with z. 
Such a collimation regime is designated in ref. |7] as the equilibrium regime. 

The general solution at k < 2 also expands as Z'^/'^ but very long wave oscillations are su- 
perimposed on this expansion [7J. These are in fact free oscillations around the equilibrium state 
arising if the fiow was injected not in the equilibrium state. They do not affect the stability of the 
fiow. The amplitude of these oscillations could be found by matching to the near zone solution at 
Z ~ 1. The spatial period of these oscillations increases with the distance as Z'^/'^. 
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The Lorentz factor of the flow is found from Eq. (H9l) . For the smooth expansion described by 
Eq. (|5T|) one can neglect the first term in the left-hand side, which yields 

2/3" 



77in 

RRl 



(52) 



The solution to this equation could be presented as 



7 = a{^)R 



(53) 



where the function ci(\E') satisfies the equation 



1 + 



2 \ 2/3 



2 

a . 



(54) 



In relativistically hot, nonmagnetized flows, the Lorentz factor grows proportionally to the flow 
radius 

7 = 7in^/i?in; (55) 

which is followed directly from the Bernoulli equation. The last expression is also valid well inside 
the light cylinder of the Poynting dominated flow because at i? <C 1, the magnetic field is nearly 
force-free and the rotation velocity is small so that one can neglect the centrifugal acceleration. 
In this case the plasma just flows along the expanding field lines and the expression ( 155|) ensures 
that the flow becomes highly relativistic, 7^1, already at i? = 1 provided the radius of the 
central body is much less than the light cylinder radius. This implies 7in ^ Rin] then Eq. 
yields a = jin/Rm- Now one sees that the hydrodynamical expression f l55|) is valid also for the 
magnetized flow in case k < 2. This in fact means that the thermal energy is transformed into the 
kinetic energy of the flow. 

It follows from Eqs. f l5B]) . (EI]) and flSSj) that the fluid could be considered as relativistically 
hot only at distances less than 

1/k 



3 '° 



-Rin7ii 



4/3 



(56) 



Beyond this distance the shape of the flow is still described by Eq. (15 ip whereas the Lorentz factor 
grows as 7 = i? until the flow ceases to be Poynting dominated [7] . 



4.3 The case k, = 2. 

If K = 2, the solution to Eq. fHTl) is written as [3l [7] 

3 



R 



/3-1/4 



1/4 



^ > 1/4- 



(57) 
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Let us consider these two cases separately. 

The /3 > 1/4 solution is similar to the equilibrium solution ( ISTl) and goes to this solution as 
f3 ^ 1/4. Note that the solution ( l5Ti) goes to F oc Z^/^ as k — i- 2 therefore the solution ( 1571) could 
be considered limiting case k — )■ 2. 

Substituting the solution ( 157|) to Eq. (149|) . one gets the equation for the Lorentz factor of the 
flow. Just as in the previous section, the solution could be presented as 

7 = ay/Z; /3 > 1/4; (59) 

where a satisfies the equation 



i + -r 



2/3 

° ^ (60) 



^ v/3(/3 - 1/4) 



The condition 7in ^ -Rm (see the previous section) implies 



3^V?EV4 7in .... 



So in this case the acceleration regime is close to that for the equilibrium fiow ( 155|) extrapolated 
to K = 2. 

Now let us consider the /3 < 1/4 solution fISH]) . The constant C in this solution is not defined; it 
could be found only by matching to the near zone solution. If the fiow was not coUimated at Z ~ 1, 
there should be C ~ 1. The /3 < 1/4 fiow is coUimated slower than Y oc Z^/^; this solution could 
be obtained by neglecting the second term in the left-hand side of the governing equation (H7|) and, 
correspondingly, in the tramsfield equation f l25p . Such a fiow was called in ref . [7] non-equilibrium 
because in this case, the pressure of the poloidal magnetic field is negligible small so that the fiow 
could be conceived as composed from coaxial magnetic loops. 

In order to find the Lorentz factor in the case /3 < 1/4 let us substitute the solution ( 158|) to 
Eq. f HU]) . One sees that one can neglect the second term in the left-hand side. Inspection of the 
obtained equation shows that 7 should grow slower than R so that one can eventually neglect the 
second term in the square brackets. Then one finally finds 



1 



This asymptotics coincides with that for the cold fiow [31 [7], which in fact means that the thermal 
energy is not transformed into the kinetic energy. In the relativistically hot fiow, the specific 
plasma energy fiux could be presented, with the aid of Eq. ( 135|) . as /i7 oc (7/r)^/^. In the non- 
magnetized plasma, /i7 = const therefore 7 oc r. In the case under consideration, /i7 decreases 
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with the distance, which means, according to the energy equation (jTTj) that the plasma thermal 
energy is converted into the Poynting flux. 

According to Eqs. (I36l) . fl58|) and fl62|) . the flow becomes cool at the distance 

-| 2/(3+VT=4^) 

(63) 



'(3 /ifnTin 



C 

Beyond this point, the flow shape and the Lorentz factor are still described by Eqs. fl58|) and (|62|) . 
correspondingly until the flow ceases to be Poynting dominated [7] . 



4.4 The case k> 2. 

At K > 2, the solution to the governing equation (H7|) goes, at large Z, to a linear function 

R = eZ; (64) 

which means that the flow becomes radial at large distances. The flow could be collimated, R Z, 
if K only slightly exceeds 2 or/and (3 is large; the final coUimation angle could be presented as[7] 
(note that our definition of G differs from that in ref. [?] by a factor 3^^^) 



e 



31/4 



K 



1/[2(k-2)] 



(65) 



'71 \K-2J V (3 

where T{x) is the gamma function. One sees that the coUimation angle rapidly increases with 
increasing k and decreasing /3: 6 = 0.013//32-5 at k = 2.2, 9 = 0.26//3 at k = 2.5 and 9 = 0.74/v^ 
at K = 3. Thus the flow could be collimated only at /« < 3. 

One can get a simpler expression for the coUimation angle assuming that k — 2 is a small number 
and making use of Stirling's approximation for the gamma function; the resulting expression reads 
as 

121/4 

9 = 777-^. (66) 



One can check that even at k = 3 this expression provides a 10% accuracy. 
It has been shown in ref. [7j that the coUimation occurs only in the region 



Z<Zo 



K-2 



(67) 



In this region, the shape of the jet is described by Eq. fl5T|) whereas the Lorentz factor grows 
according to Eq. fl55|) . Therefore the k < 2 scalings remain valid at k > 2 but only till some 
limiting distance. At Z > Zq, the jet shape approaches the cone flM|) . The curvature of the flux 
surfaces decreases in this region; it could be estimated from the governing equation ( H7I) as 

cPR f3Q 



(PZ 



(6^ 
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Figure 1: Sketches of the jet radius (a) and Lorentz factor (b) as functions of the distance from 
the origin. The curve 1 describes the case k < 2; both the jet radius and the Lorentz factor grow 
as a power law with the index less than 1/2, see Eqs. (1511) and fl55|) respectively. The case k = 2, 
/3 > 1/4 may be considered as the limiting case of this regime with both the radius and the Lorentz 
factor growing oc ^J~Z. The case k = 2, /3<l/4is shown by curves 2; the jet expands according 
to Eq. (!58|) faster than ^J~Z whereas the Lorentz factor grows according to Eq. (162!) slower than 
a/Z. The curve 3 is for k slightly larger than 3; at the initial stage (in the region ( 167|) ) the jet 
is rapidly expanded (according to Eq. (EH)) ai^^ accelerated (according to Eq. ( |55l) ). Then it 
becomes conical preserving the acquired collimation angle (in logarithmic coordinates, this looks 



Then Eq. (I49l) for the Lorentz factor is written as 



*0 

i + -r 



77in 

QZRl 



2/3' 



K-2 



+ 



(69) 



The second term in the right-hand side rapidly decreases with Z and becomes less than the first 
term at 



Z>Z^= (/304) 
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/3 



1/(k-2). 



(70) 



we used Eq. (166|) for 9 in the last equahty. Since Eq. ( I69l) is valid only at Z > Zq and the ratio 

2/(k-2) 



1 2/(4-k){k-2) 

2 



(71) 



is less than unity at 2 < k < 3, one can neglect the second term in the right-hand side of this 
equation. Moreover one can check that 7 should grow slower than R oc Z in this case so that one 
can eventually neglect the second term in the square brackets. Then one finds that the Lorentz 
factor grows as 



7 



1 



(72) 



which coincides with the corresponding scaling for cold flows (see Eq. (112) in ref.[7]). As we have 
already discussed at the end of sect. 4.3, this implies that the thermal energy is transformed into 
the Poynting flux, not into the kinetic energy. 

One should stress that the scaling fl72l) is obtained for the conical part of the jet where the 
curvature of the flow lines, d'^R/dZ'^, is determined by small deviations from the straight line. 
Therefore Eq. fl72]) is valid only if R{Z) could be found from the governing equation fHTl) with the 
necessary accuracy. We show in the Appendix that this is indeed the case. 

All the above scalings were obtained at the condition that the plasma is relativistically hot, 
h ^ 1. However, the expanding flow continuously cools down so that eventually this condition is 
violated. In the region Z < Zi, the Lorentz factor increases as 7 = 7in-R/-Rin therefore Eq. f l35p 
implies that the flow remains hot till Z ^ Z\ provided 



kx-1/{4-k) 



(73) 



Let us assume that this condition is fulfilled; then the transition to the cold flow occurs in the 
region Z > Zx where the Lorentz factor grows according to Eq. fl72l) . It follows from Eqs. fl36|) . 
(1^ and (1721) that the transition to the cold flow occurs at the distance 

2/(k+2) 



(74) 
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The flow Lorentz factor at the transition point is 

7==(^)"'""'(«?»7,>y'"-^'""*^ (75) 
With account of the condition (173!) . one gets 

7c ^ f R -m 

Irnhn V 7in 

Taking into account that i?in < 7in in relativistically hot flows, one sees that 7c < 7in^in, which 
means that most of the thermal energy is transferred not to the kinetic energy but to the Poynting 
flux. 

When the plasma becomes cool, the flow still remains Poynting dominated. It has been shown 
in ref. [7] that at k > 2, the cold flow could be accelerated according to the scaling (!72|) only till 
a terminal Lorentz factor 7f ~ (7max/0^) after which the acceleration becomes extremely slow, 
7 ~ 7( {hiZ / Zt)^^^ . Therefore efficient conversion of the Poynting flux into the flow kinetic energy 
could occur only if 7max ^ 7t, which is equivalent to the condition 

7maxe < 1. (77) 

Now we see that this statement could be generalized to hot flows. Since the cold flow could not be 
accelerated efficiently to 7 larger than 74, the final Lorentz factor is the smallest of 7c, 7* and 7max- 
If the flow was initially Poynting dominated, /iin7in ^ 7max, it becomes cool still being Poynting 
dominated, 7c ^ 7max? and therefore the efficient conversion of the Poynting flux into the kinetic 
energy could anyway occur only if the condition (1771) is fulflUed. 




5 Conclusions 

In this paper, we studied hot, strongly magnetized, relativistic jets at large distances from the 
source. The acceleration zone of relativistic, Poynting dominated jets spans a very large range of 
scales therefore the processes far away from the source are of special interest. Multi-scale systems 
generally pose a strong challenge to numerical simulations. On the other hand, they are suitable for 
asymptotic analysis. We derived equations governing the flow in the far zone; they do not contain 
terms that nearly cancel each other, as the general MHD equations do, and therefore they could 
be solved relatively easily. We concentrated on relativistically hot flows because cold, Poynting 
dominated jets has already been extensively studied. Our results directly generalize the results of 
ref. [7] for cold flows. 

Making use of the obtained asymptotic equations, we studied in detail the structure of the 
jet from a rigidly rotating source, Vt = const. The coUimation and acceleration of the flow are 
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intimately connected and determined by the distribution of the confining pressure. We adopted a 
power law distribution, pext oc z~'^, so that the solution depends on two parameters, the index k 
showing how fast the pressure decreases, and the parameter /3 defined by Eq. fj48|) and showing 
how strong is the external pressure extrapolated to the light cylinder. We have found that hot, 
Poynting dominated jets are collimated exactly as in the cold case. Namely, at k < 2 the fiow 
shape is described by a power law function ( Eq. ( 15T]) for k, < 2 and Eqs. ( [57|) and ( 158]) for /t = 2) 
so that the jet opening angle continuously decreases. At k > 2 the fiow becomes asymptotically 
radial, the final collimation angle being presented by Eq. ( !66l) . 

When the fiow expands, the fiuid cools down. We have shown that at k < 2 or k = 2, 
P > 1/4, the thermal energy is converted into the kinetic energy of the fiow. On the contrary, at 
K = 2, P < 1/4 or K > 2, the thermal energy is converted into the Poynting fiux so that even the 
relativistically hot fiow is accelerated as if it is cold. Note that acceleration regimes of cold fiows 
are also different in these parameter ranges. At k < 2, the cold fiows are accelerated until the 
equipartition between the Poynting and the plasma kinetic energy fluxes is eventually achieved. 
At K > 2, the flow acceleration ceases when a limiting Lorentz factor is achieved so that the flow 
could remain Poynting dominated [7]. One now sees that this conclusion remains valid also for 
relativistically hot but Poynting dominated flows. 

In gamma-ray bursts, the relativistic jet is formed during the collapse of star's core. In this 
case, the outflow is initially relativisticallly hot. Observations of the burst afterglows suggest that 
the flnal opening angle of the jet is a few degrees whereas the flnal Lorentz factor is at least a few 
hundreds. This implies O7 ^ 1, which is characteristic for the case k, > 2. We have shown that 
in this case the thermal energy is converted into the Poynting flux, not into the kinetic energy. 
Therefore the thermal acceleration could not help to transform the Ploynting flux into the plasma 
energy. Some sort of magnetic dissipation is necessary in order to utilize the electro-magnetic 
energy completely. 

This work was supported by the US-Israeli Binational Science Foundation under grant number 
2006170 and by the Israeli Science Foundation under grant number 737/07. 

Appendix. Corrections of the order of hj/fi to the flux 
surface shape. 

We found the jet structure neglecting hj as compared with /i in the Bernoulli equation fl26|) . 
Then the shape of the flux surfaces is given by Eq. fH5]) where R{Z) satisfles the governing 
equation f H7|) . The accuracy of the procedure, to within a factor h'-f/fi, is generally sufficient in 
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the Poynting dominated jets. However, a special attention should be given to the case k> 2 when 
the flow becomes radial at large distances (see sect. 4.4). The Lorentz factor of the flow depends, 
according to Eq. f l49|) . on the curvature of the flow line, (PR/dZ"^. When the flow is close to radial, 
the curvature is determined by small corrections to the function f lM|) . In sect. 2.4, we found the 
Lorentz factor of the flow in the form ( 172|) assuming that the curvature of the flow lines could 
be found from the solution to the governing equation. Inasmuch as this equation is valid only to 
within the factor /i7//i, we have to check whether the curvature, and consequently the Lorentz 
factor, thus found is correct. To do so we substitute the Lorentz factor from Eq. (I49l) into the 
Bernoulli equation (l26ll and find the corresponding corrections to the shape of the flux surfaces. 
The expression f l72|) is valid while the curvature due to this corrections remains small as compared 
with the curvature fl68l) obtained from the solution of the governing equation. 
Let us present the shape of the flux surfaces as (cf. Eq. fH3l) ) 



X(Z, = RiZ)^Jl + 6{Z, M/)]; (78) 

where 5(\1/, Z) describes corrections to the shape of the flux surfaces due to a non-zero h'y/fi. 
Substituting this expression into Eq. fl26|) and linearizing with respect to small 6 and hj/fi, one 
gets, with the aid of Eq. fl39|) . 

— = • (79) 

where 

S = ^/^o; 7max = /i(^o) = 2n^<i!o/v- (80) 
With 7 from Eq. fl72|) and h from Eq. fl35l) . one finds 

dS 2e4/3/31/3^^^^57/3^{4-K)/3' ^ ) 

which yields 

o 7 1/3 p2/3 

r 'J"-in7in -^in /qoN 

8eV3/3V3^^^^^4/3^(4-.)/3- ^^^) 

Substituting this expression into Eq. ( ITHj) and differentiating twice with respect to Z, one finds 
the curvature of the flux surface as 



1/3 d2/3 



d^X _ ^ d^R (4-/t)(7-/t)/zin7i;^X 

dZ^ [dZ^ 2401/3/31/3^^^^^4/3^(7~K)/3 



(83) 



Here we take into account that R ~ QZ. 

The Lorentz factor of the flow could be determined from the solution to the governing equation 
only if the second term in brackets in Eq. fl83|) is small as compared with the first one. The 
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governing equation yields Eq. (|68l) for cpR/dZ"^ therefore the ratio of the first to the second terms 
is written as 

24(/3e)^/37^,,zf(5-2-) 



(4-«:)(7-«:)/^in7ifi?f 



(84) 

At the point Zq beyond which the acceleration law (1721) is valid, this ratio could be estimated, 
with the aid of Eqs. (166|) and ( 1671) . as 



210/334/3 



(4- k)(7- k) 



1 



e \ K 



5-2.1 4/3(«-2) 



7max ^gg^ 



, 1/3 p2/3 ■ 
"■inTin -'^in 



Taking into account that /iin7in ^ 7max and 7in ^ -Rin one sees that this quantity is large. If 
K < 2.5, the ratio (18^ remains large aX Z > Zi therefore in this case the governing equation 
correctly describes the curvature of the filed line, which justifies the use of Eq. (!72|) for the 
Lorentz factor of the flow. At k > 2.5, the above ratio decreases with the distance therefore 
eventually Eq. (1721) could become invalid; we do not consider this case because we are interested 
in strongly collimated jets whereas according to Eq. fl66|) . the flow is only weakly collimated at 
K > 2.5. 

An important point is that the above conclusions were obtained at the condition that the flow 
is relativistically hot. It is shown in ref. [7] that in cold flows, the corrections of the order of 7//^ 
to the flux surface shape eventually become significant in /t > 2 flows; then the acceleration stops. 
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